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Abstract 

A status report on the theory and phenomenology of rare radiative B decays in the 
standard model is presented with emphasis on the measured decays B — ► X s ^ and 
B — > K*j. Standard model is in agreement with experiments though this comparison 
is not completely quantitative due to imprecise data and lack of the complete next-to- 
leading order contributions to the decay rates. Despite this, it is possible to extract non- 
perturbative parameters from the shape of the photon energy spectrum in B — > X s ^, 
such as the 6-quark mass and the kinetic energy of the b quark in B hadron. The 
measured decay rate B(B —y X s + 7) = (2.32 ± 0.67) x 10~ 4 can also be used to extract 
the CKM ratio |T4 s |/|V cfc |, yielding | V^ b | = 1.10 ± 0.43. Issues bearing on the 
determination of the parameters of the CKM matrix from the CKM-suppressed decays 
B — ► Xd + 7 and B — > (p, uj) + 7 are also discussed. It is argued that valuable and 
independent constraints on the CKM matrix can be obtained from the measurements 
of these decays, in particular those involving neutral £>-mesons. 
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1 Introduction 



The last two years have seen the first observations of the electromagnetic penguins in B 
decays by the CLEO collaboration. These include the measurements of the exclusive decay 
rate, B(B -> K* + 7) = (4.5 ± 1.0 ± 0.9) x 10~ 5 [|, and the inclusive rate B(B -> X s + 7) = 
(2.32 ±0.67) x 10~ 4 §, yielding R(K*/X S ) = T(B -> K* + i)/T(B -> X s ± 7 ) = 0.19 ±0.09. 
In addition, the charged and neutral 5-meson decay rates are found equal within experimental 
measurements. The inclusive decay rate is in agreement with the predictions of the standard 
model H|, H, with rate estimated as B(B -> X, + 7) = (2.55 ± 1.28) x 10~ 4 assuming 
|Vts|/|K;&| = 1.0 M. Conversely, one can vary the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
element ratio and determine it from BiB — > X s ± 7), which yields |V^ S |/|V^,| = 1.10 ± 0.43. 
The ratio R(K*/X S ) is well explained by the QCD-sum-rule based estimates of the recent 
vintage |7], || and by wave-function models combined with vector meson dominance (local 
parton-hadron duality) 0. This and the near equality of the charged and neutral decay rates 
imply that the observed radiative B decays are dominated by the (common) electromagnetic 
penguin (short distance) amplitudes and the contributions from I?-meson-specific diagrams 
(weak annihilation, IV ± -exchange) are small. 

The photon energy spectrum in B — > X s ± 7 yields information on the structure function 
of the photon in the electromagnetic penguins ||, [10], IT], [T2]] . In specific models I3J , 



this information can be transcribed in terms of non-perturbative parameters, such as the 
6-quark mass and the kinetic energy of the b quark in B hadron. These quantities can also 
be estimated in the framework of the heavy quark effective theory |15j combined with QCD 



sum rules [0, [171 ■ The results of a recent analysis of the CLEO data on B — > X s + 7 || 
are consistent with the values expected from such theoretical considerations. However, this 
agreement is presently not completely quantitative due to imprecise data. 

There is considerable interest in measuring the CKM-suppressed radiative B decays, such 



as B — > Xd ± 7 |L8| and B — > (p, u) + 7 0. A determination of the CKM parameters from 
eventual measurements of these decays requires careful treatment of the competing short- 
distance (SD) and long-distance (LD) effects. This problem can be formulated in terms 
of model-independent correlation functions involving matrix elements of a few dimension-6 
operators in an effective theory. Techniques, such as the QCD sum rules, can then be invoked 
to estimate them. In (TP], ^TJ, the leading LD-effects in the exclusive decays B — > (p, uj) + 7 
are calculated in terms of the weak annihilation amplitudes. The largest such effects may 
show themselves in the charged S ± -decays, B ± — > + 7, contributing up to 0(15%) of the 
corresponding SD-amplitudes; their influence in the neutral 5-decays is estimated to be much 
smaller. Hence, there are good theoretical reasons to plead that the decays B° —> (p°, uj) + 7 
and B —>■ Xd ± 7 are well suited to determine the CKM parameters. We take up these issues 
in this status report. 



2 Estimates of B(B -> X s + 7) in the SM 

The framework that is used generally to discuss the decays B — » X s ± 7 is that of an effective 
theory with five quarks, obtained by integrating out the heavier degrees of freedom, which 
in the standard model are the top quark and the H^-boson. A complete set of dimension-6 
operators relevant for the processes b — > s + 7 and b ^ s + j + g is contained in the effective 
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where Gp is the Fermi constant coupling constant, Cj(fi) are the Wilson coefficients evalu- 
ated at the scale /x, and Xt = Vt b V t * with Vij being the CKM matrix elements. The overall 
multiplicative factor At follows from the CKM unitarity and neglecting X u . The operators Oj 
read 
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where e and g s are the electromagnetic and the strong coupling constants, respectively. 
In the magnetic moment type operators 7 and 8 , F^ u and G^ u denote the electromag- 
netic and the gluonic field strength tensors, respectively. The subscripts on the quark fields 
L = (1 — 7s)/2 and R = (1 + 7s)/2 denote the left and right-handed projection operators, 
respectively. QCD corrections to the decay rate for b — > 57 bring in large logarithms of the 
form a™(mw) ^og m (m b /M), where M = m t or mw and m < n (with n = 0, 1, 2, ...). Using the 
renormalization group equations the Wilson coefficient can be calculated at the scale \l ~ m b 
which is the relevant scale for B decays. To leading logarithmic precision, it is sufficient to 
know the leading order anomalous dimension matrix and the matching Ci(fi = mw) to lowest 
order (i.e., without QCD corrections) [plfl . The 8x8 anomalous dimension matrix is given in 
from where references to earlier calculations can also be obtained, the Wilson coefficients 
are explicitly listed in || and the numerical values of these coefficients being used here can 
be seen in ||. 

It has become customary to calculate the branching ratio for the radiative decay B — > 
X s + 7 in terms of the semileptonic decay branching ratio BiB — > Xii/f) 
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where, in the approximation of including the leading-order QCD correction, T s i is given by 
the expression 
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The phase space function g(z) and the function f(z) due to one-loop QCD corrections can 
be seen in [[□[]. The radiative decay rate T(B — > X s + 7) are worked out in [|3], [J, taking into 
account 0(a s ) virtual and bremsstrahlung corrections. In calculating the matrix elements in 
these papers, the on-shell subtraction prescription for the quark masses has been used. Due 



Parameter 


Range 


rn t (GeV) 
a (GeV) 
A 5 (GeV) 
B(B -> 
m c /m b 
m w (GeV) 

a QED 


170 ±11 
5.0+1° 
0.195i°';°f 
(10.4 ±0.4)% 
0.29 ±0.02 
80.33 
130.0 



Table 1: Values of the parameters used in estimating the branching ratio BiB — » X s + 7) in 
the standard model. 

to the explicit factors of the running quark masses in the operators O7 and 8 , the m^-factor 
contained in the decay rate F(B — > X s + 7) is replaced by the following product 

ml — > m b (pole) 3 m b (fj,) 2 , (5) 

where m b (pole) and m fe (/i) denote the pole mass and the MS-running mass of the b quark, 
respectively. Since, in the leading order in a s , the semileptonic decay width T s i depends 
on the product mb(pole) 5 , the ratio of the two decay widths brings in the correction factor 
R(m b ,fi): 

R(m b ,fi) = [m b ({i)/m b (pole)} 2 , (6) 
as also remarked in p|. At the one-loop level, these masses are related: 

m b{n) _ r as(fj) 14/fth 4 a 8 (/A 
m b (pole) a s (m b ) 3 i 

where /?o = 23/3. The parameters used in estimating the inclusive rates for B(B — >• X s + 7) 
are summarized in table 0. 

We now discuss BiB — > X s + 7) in the standard model and theoretical uncertainties on 
this quantity 0. 

• Scale dependence of the Wilson coefficients. 

The largest theoretical uncertainty stems from the scale dependence of the Wilson coefficients. 
As derived explicitly in |J, the decay rate for B — > X s +j depends on seven of the eight Wilson 
coefficients in H e ff(b —> s), once one takes into account the bremsstrahlung corrections and 
is not factored in terms of a single (effective) coefficient, namely C7 , that one encounters 
for the two-body decays 6 s + 7 |, |. Numerical values of the two dominant effective 
coefficients, Cy" and , as one varies /i, the QCD scale A 5 , and the (running) top quark 
mass in the MS-scheme fht(mt) in the range given in table 1, are: 

C 7 e// = C 7 - - C 6 = -0.306 ± 0.050, 

Ct ff = C 8 + C 5 = -0.146 ±0.020. (8) 
This is the dominant theoretical error on BiB — ► X s + 7), contributing about ±35%. 

• Scale-dependence of m b (fi) in the operators O7 and 0%. 
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This brings into fore the extra (scale-dependent) multiplicative factor -R(mj,, /i) for the branch- 
ing ratio BiB —>■ X s + 7), as discussed above. Intrinsic uncertainties in the concept of the 
pole mass due to infrared renormalons suggest that one should express all physical results in 
terms of the running masses p2| . This requires recalculating the decay rate BiB — > X s + 7) 
with the running masses, incorporating resummations of the kind recently undertaken for the 
semileptonic B decay rates |23| . 



• Extrinsic errors in BiB — > X s + 7) 

The next largest error arises from the parameters which are extrinsic to the decay B — > X s + 7 
and have crept in due to normalizing the branching ratio BiB — > X s + 7) in terms of BiB — > 
Xlvi). The first of these extrinsic errors is related to the uncertainty in the ratio m c /m{,. 
Using for the b quark pole mass m^pole) = 4.8 ± 0.15 GeV |16| and m^, — m c = 3.40 GeV 
f24j| , one gets m c /mb = 0.29 ± 0.02. Taking into account the experimental error of ±4.1% on 
BiB — > Xtvi) pij] , one estimates an extrinsic error of ±12% on BiB — ► X s + 7). 

Assuming |T4s|/|K:& = 1 HHj the branching ratio BiB — > X s j) calculated as a function of 
the top quark mass is shown in Fig. 1 ||. For all three solid curves the quark mass ratio is 
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Figure 1: BiB — > Xg'y) as a function of the MS top quark mass. The three solid lines corre- 
spond to the variation of the parameters fi and A5 as described in the text. The experimental 
(±l<r)-bounds from CLEO [0] are shown by the dashed lines (from 0). 

fixed at m c /mb = 0.29. The top solid curve is drawn for /i = 2.5 GeV and A 5 = 0.260 GeV. 
The bottom solid curve is for /i = 10 GeV and A5 = 0.145 GeV, and the middle solid curve 
corresponds to the central values of the input parameters in table 1. Using m t = (170 ± 11) 
GeV, and adding the extrinsic error, one obtains: 



BiB -> X s + 7) = (2.55 ± 1.28) x 10" 



(9) 



to be compared with the CLEO measurement BiB — > X s ± 7) = (2.32 ± 0.67) x 10~ 4 . The 
(±l<r)-upper and -lower bound from the CLEO measurement are shown in Fig. 1 by dashed 
lines. The agreement between SM and experiment is good, given the large uncertainties on 
both. In 0, the branching ratio BiB — » X s ± 7) has been calculated as a function of the 
CKM matrix element ratio squared (| V ts \/\ Kfe|) 2 , varying m t , H and A 5 in the range specified 
in table 1. Using the (±lcr)-experimental bounds on BiB — > X s ± 7), one infers |J: 



Vj/Wr 



cb\ 



1.10 ±0.43 



(10) 
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which is consistent with the indirect constraints from the CKM unitarity [^6] yielding | V^|/| V^,| : 
1.0 but imprecise. Further improvements require reducing the perturbative scale (/^-dependence 
of the decay rate, which in turn implies calculations of the next-to-leading order terms, and 
more accurate measurements. 



3 Photon Energy Spectrum in B —> X s + 7 

The two-body partonic process b — > sj yields a photon energy spectrum l/TdT(b — > 57) = 
8(1— x) , where the scaled photon energy x is defined as E 1 = (m 2 — m 2 ) / (2 m&) x; x then varies 
in the interval [0, 1]. Perturbative QCD corrections, such as b — > 57 + g, give a characteristic 
bremsstrahlung spectrum peaking near the end-points, J5 7 — > E™ ax and E 1 — > 0, arising 
from the soft-gluon and soft-photon configurations, respectively. As long as the s-quark 
mass is non-zero, there is no collinear singularity in the spectrum. Near the end-points, one 
has to improve the spectrum obtained in fixed order perturbation theory. This is usually 
done by isolating and exponentiating the leading behaviour in a em a s (/i) m log n (l — x) and 
a em a s (ii) m log n x, with m < n, where /x is a typical momentum in the decay B — > X s + 7. 
The running of a s is a non-leading effect, but as it is characteristic of QCD it modifies the 
Sudakov-improved end-point photon energy spectrum |ll] compared to its analogue in QED 
27f . Away from the end-points, the photon energy spectrum has to be calculated completely 



in a given order in a s in perturbation theory || ||. 

The complete photon energy spectrum in B — > X s + 7 is at present not calculable in QCD 
from first principles. The situation is very much analogous to that of other hadronic structure 
functions. It has been observed in a number of papers |], [II], [TTJ , that the x-moments of the 



inclusive photon energy spectrum in B — > X s + 7 and those of the lepton energy spectrum in 
the decay B — > X u lvi are related. Defining the moments as: 

I f M B /m b dY 

M n {B^X s + 1 ) = -J dxx n ~ l — (11) 

pM B /m b d 1 dYt 

M n (B^XJu e ) = - dxx n —{--±) 

Jo ax I £ ax 

n r M B/m b n ,dr e 
dxx 



T e Jo dx 

the ratios of the moments are free of non-perturbative complications. The moments M. n have 
been worked out in the leading non-trivial order in perturbation theory and the results can 
be expressed as: 

M n ~ 1 + — C F (A log 2 n + B logn + const.) (12) 
2ir 

where Cf = 4/3, the leading coefficient is universal with A = — 1 ||27|| , and the non-leading 
coefficients are process dependent; B = 7/2 | and B = 31/6 |2l| , for B — > X s + 7 and 
B — > X u £ui, respectively. Measurements of the moments could eventually be used to relate 
the CKM matrix element V ts and V u b. That this method will give competitive values for V u b, 
however, depends on whether or not the coefficient functions in T(B — >• X s + 7) discussed in 
the previous section are known to the desired level of theoretical accuracy. 

We shall leave such theoretically improved comparisons for future Rencontres de Moriond 
and confine ourselves to the discussion of the present state-of-the-art comparison of the mea- 
sured photon energy spectrum in B — > A s +7 with the perturbative QCD-improved treatment 
of the same. The analysis that we discuss here B B[ treats the non-perturbative effects in 



terms of a 5-meson wave function. In this model |13j , which admittedly is simplistic but not 
necessarily wrong, the b quark in B hadron is assumed to have a Gaussian distributed Fermi 
motion determined by a non-perturbative parameter, pp, 

4 — p 2 

4>\P) = r- 3 eX P( 2") ' P=\P\ ( 13 ) 

V^Pf* Pf 2 

with the wave function normalization / °° dpp 2 <p(p) = 1. The photon energy spectrum from 
the decay of the I?-meson at rest is then given by 

dr rVmax „ (TV, 

wrl * r * ) 4 KfA) • (14) 

where p ma x is the maximally allowed value of p and j^r- is the photon energy spectrum from 
the decay of the 6-quark in flight, having a momentum-dependent mass W(p). 

An analysis of the CLEO photon energy spectrum has been undertaken in || to determine 
the non-perturbative parameters of this model, namely Tadpole) and pp. The experimental 
errors are still large and the fits result in relatively small x 2 values; the minimum, Xmin = 
0.038, is obtained for p F = 450 MeV and m b (pole) = 4.77 GeV, in good agreement with 
theoretical estimates of the same, namely m^pole) = 4.8 ± 0.15 GeV Jl6[ and pp = fi 2 /2 = 



0.25 ± 0.05 GeV 2 obtained from the QCD sum rules [17|. In Fig. 2 we have plotted the 
photon energy spectrum normalized to unit area in the interval between 1.95 GeV and 2.95 
GeV for the parameters which correspond to the minimum x 2 (solid curve) and for another 
set of parameters that lies near the x 2 -boundary defined by x 2 — Xmin + 1- (dashed curve). 
Data from CLEO @ are also shown. Further details of this analysis can be seen in 0. 




E 7 [GeV] Fig . 4a 



Figure 2: Comparison of the normalized photon energy distribution using the CLEO data 
corrected for detector effects and theoretical distributions from , both normalized to 
unit area in the photon energy interval between 1.95 GeV and 2.95 GeV. The solid curve 
corresponds to the values with the minimum % 2 , (m 9 ,pi?) = (0,450 MeV), and the dashed 
curve to the values (m g ,p F ) = (300 MeV, 310 MeV). 



4 Inclusive radiative decays B —> + 7 

The theoretical interest in the standard model in the (CKM-suppressed) inclusive radiative 
decays B —>■ X^+'y lies in the first place in the possibility of determining the CKM-Wolfenstein 



parameters p and 77 [30]. The relevant region in the decays B — > X^+7 is the end-point photon 
energy spectrum, which has to be measured requiring that the hadronic system X^ recoiling 
against the photon does not contain strange hadrons to suppress the large- E 1 photons from 
the decay B — > X s + 7. Assuming that this is feasible, one can determine from the ratio of 
the decay rates BiB — > Xd + r y)/B(B — > X s + 7) the CKM-Wolfenstein parameters. This 
measurement was proposed in [fjjjfl , where the final-state spectra were also worked out. 

In close analogy with the B — > X s +7 case discussed earlier, the complete set of dimension- 
6 operators relevant for the processes b — > and b — > djg can be written as: 

AC 8 

H Q, _> d ) = —i-it EQW^). (15) 
V2 J=1 

where = V^VL for j = t,c,u. The operators O,-, j = 1,2, have implicit in them CKM 



factors. In the Wolfenstein parametrization |p0| , one can express these factors as : £ u = 
AX 3 (p — ir]), £ c = — v4A 3 , & = — — £ c . We note that all three CKM-angle-dependent 
quantities £j are of the same order of magnitude, 0(A 3 ), where A = sin#c* ~ 0.22. This 
is an important difference as compared to the effective Hamiltonian 7i e ff(b — > s) written 
earlier, in which case the effective Hamiltonian factorizes into an overall CKM factor At. For 
calculational ease, this difference can be implemented by defining the operators 0\ and O2 



entering in H e ff(b — > d) as follows [0| : 



01 = - — {c L p^b La )(d La ^^c L p) - —(u L/ 3'j ,1 b La )(d La 'y^u Lf 3), 

2 = -7 £ (cL a 7 / "6 ia )((iL /3 7 At CL /3 ) - -^(u La Yb La )(d L ^^u L p), (16) 

and the rest of the operators (Oj; j = 3. ..8) are defined like their counterparts Oj in H e ff(b — > 
s), with the obvious replacement s — > d. With this definition, the matching conditions 
Cj(mvy) and the solutions of the RG equations yielding Cj(p) become identical for the two 
operator basis Oj and Oj. It has been explicitly checked in the 0(a s ) calculations of the 
decay rate and photon energy spectrum involving b — > dg and b — > ^7 transitions that the 
limit m u — > for the decay rate T(B — > X^ + 7) exists [0. From this it follows that, in 
the leading order QCD corrections, there are no logarithms of the type a s log(m^/m^) |2"P| . 



Some papers, estimating LD-contributions in radiative B decays, seem to contradict this by 
assuming light-quark contributions which have such spurious log- dependence. There is no 
calculational basis for this assumption. In higher orders, such terms must be absorbed in the 
non-perturbative functions. On the other hand, as far as the dependence of the decay rate 
and spectra on the external light quark masses is concerned, one encounters logarithms of the 
type a em a s (l + (1 - x) 2 ) / x\og(m 2 b /m 2 s ) (for b -> sg^) and a em a s (l + (1 - x) 2 )/x \og(ml / m 2 d ) 
(for b — > dg^f) near the soft-photon (x — > 0) region |J, which can, however, be exponentiated 
||12|| . The essential difference between T(B — > X s + 7) and T(B — > X^ + 7) lies in the 
matrix elements of the first two operators Oi and O2 (in H e ff(b — > s)) and Oi and O2 (in 
H e ff(b —> d)). The derivation of the inclusive decay rate and the final-state distributions in 
B —>■ Xa + 7 otherwise goes along very similar lines as for the decays B —>■ X s + 7 . The 
branching ratio B(B — > X^ + 7) in the SM can be written as: 

S( J B^X d + 7 ) = J D 1 |6| 2 

1 — P V D4 

{1 ~ (1 - p) 2 + r/ 2 ° 2 ~ (l-p) 2 + r] 2 ° 3+ [l-pf + r]^ 1 (17) 
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where the functions Di depend on the parameters listed in table 1. The uncertainty on this 
branching ratio from the parametric dependence is very similar to the one worked out for 
B(B — > X s + 7). For the central values of the parameters in table 1, one gets : D\ = 
0.21, D 2 = 0.17, D 3 = 0.03, D 4 = 0.10. To get the inclusive branching ratio the CKM 
parameters p and 77 have to be constrained from the unitarity fits. Taking the parameters 
from a recent fit, one gets 5.0 x 10~ 3 < \£ t \ < 1.4 x 10 -2 (at 95% C.L.) f3l], yielding an order 
of magnitude uncertainty in BiB — > Xd + 7) - hence the interest in measuring it. Taking the 
central values of the fit parameters A = 0.8, A = 0.2205, 77 = 0.34 and p = —0.07 [3TJ, one 
gets B(B — > Xd + 7) = (1.7 ± 0.85) x 10 -5 , which is approximately a factor 10 -20 smaller 
than the CKM-allowed branching ratio B{B — ► X s + 7), measured by CLEO [|TJ. 



5 Estimates of B(B —> V + 7) and Constraints on the 
CKM Parameters p and 77 

Exclusive radiative B decays B — > V + 7, with V = K*,p,u, are also potentially very 
interesting from the point of view of determining the CKM parameters [p]]. The extraction 
of these parameters would, however, involve a trustworthy estimate of the SD- and LD- 
contributions in the decay amplitudes. 

The SD-contribution in the exclusive decays (B u ,Bd) — > (K*,p) + 7, B^ — > u + 7 and 
the corresponding B s decays, B s — ► (</>, K*) + 7, involve the magnetic moment operator O7 
and the related one obtained by the obvious change s — > d, O7. The transition form factors 
governing the radiative B decays B —>■ V + 7 can be generically defined as: 

(V, \\^a, u q%\B) = ie» upa ei X) PBPvF^ Y (0)- (18) 

Here V is a vector meson with the polarization vector e^, V = p,cu, K* or 0; B is a generic 
£>-meson B u , Bd or B s , and ip stands for the field of a light u, d or s quark. The vectors pb, Pv 
and q = Pb —Pv correspond to the 4- momenta of the initial 5-meson and the outgoing vector 
meson and photon, respectively. In ([18]) the QCD renormalization of the ipa^ u q v h operator 
is implied. Keeping only the SD-contribution leads to obvious relations among the exclusive 
decay rates, exemplified here by the decay rates for (B u , Bd) — > p + 7 and (B u , Bd) — > K* + 7: 

I V ts\ 

where <& U) d is a phase-space factor which in all cases is close to 1 and Ki = [Fg x ~^ pl / ' p^~* K 7 ] 2 
is the ratio of the (SD) form factors squared. The transition form factors F s are model depen- 
dent. However, their ratios, i.e. Ki, should be more reliably calculable as they depend essen- 
tially only on the SU(3)-breaking effects. If the SD-amplitudes were the only contributions, 
the measurements of the CKM-suppressed radiative decays (B u , Bd) — > p + 7, Bd — > cu + 7 and 
B s — > K* + 7 could be used in conjunction with the decays (B u , B d ) — > K* + 7 to determine 
the CKM parameters. The present experimental upper limits on the CKM ratio |Vtd|/|K s | 
from radiative B decays are indeed based on this assumption, yielding 



r((B u ,B d ) 


-»• P + 7) 






T((B u ,B d )- 


+ K* + 7) 


N 2 





\V t d\ 

t^I< 0.75, (20) 

\Vts 



with a theoretical dispersion estimated in the range 0.64 - 0.75, depending on the models used 
for the SU(3) breaking effects in the form factors |7], |33 |. 



The possibility of significant LD- contributions in radiative B decays from the light quark 
intermediate states has been raised in a number of papers [[34]]. Their amplitudes necessarily 



involve other CKM matrix elements and hence the simple factorization of the decay rates 
in terms of the CKM factors involving \Vu\ and \Vt s \ no longer holds thereby invalidating 
the relationships given above. In what follows, we argue that the CKM-analysis of charged 
.B-decays, B ± — > p ±r y, would require modifications due to the LD-contributions but the 
corresponding analysis of the neutral -B-decays B — > [p°,u)j remains essentially unchanged. 

The LD-contributions in B — > V + 7 are induced by the matrix elements of the four- 
Fermion operators 0\ and O2 (likewise 0\ and O2). Estimates of these contributions require 
non-perturbative methods. This problem has been investigated recently in |19|, |2(| using a 
technique which treats the photon emission from the light quarks in a theoretically consistent 
and model-independent way. This has been combined with the light-cone QCD sum rule 
approach to calculate both the SD and LD - - parity conserving and parity violating - 
amplitudes in the decays B u ^ — > p{oj) + 7. To illustrate this, we concentrate on the 
decays, B^ — ^ p ± + 7 and take up the neutral B decays B d — > p{oj) + 7 at the end. The 
LD-amplitude of the four-Fermion operators 0\, O2 is dominated by the contribution of the 
weak annihilation of valence quarks in the B meson. It is color-allowed for the decays of 
charged B^ mesons, as shown in fig. 3, where also the tadpole diagram is shown, which, 
however, contributes only in the presence of gluonic corrections, and hence neglected. In the 
factorization approximation, one may write the dominant contribution in the operator O2 
(here 0' 2 is the part of O2 with the CKM factor 

(fn\0' 2 \B) = (p\dr (t u\0)('y\ur»b\B) + (frf\dT^u\0) (0\uT%\B) , (21) 

and make use of the definitions of the decay constants 

(o|«r M 6|s) = i P(l f B , 

( P \dT,u\0) = e^m p f p , (22) 

to reduce the problem at hand to the calculation of simpler form factors induced by vector 
and axial-vector currents. 





Figure 3: Weak annihilation contributions in B u — > p7 involving the operators 0[ and 0' 2 
denoted by with the photon emission from a) the loop containing the b quark, b) the loop 
containing the light quark, and c) the tadpole which contributes only with additional gluonic 
corrections. 

The factorization approximation assumed in [|l^, ^] has been tested (to some extent) in 
two-body and quasi-two body non-leptonic B decays involving the transitions b —> ccs and 



b — > cud. It has not been tested experimentally in radiative B decays. From a theoretical 
point of view, non-factorizable contributions belong to either the 0(a s ) (and higher order) 
radiative corrections or to contributions of higher-twist operators to the sum rules. Their 
inclusion should not change the conclusions substantially. 

The LD-amplitude in the decay B u — > p ± + 7 can be written in terms of the form factors 
Fi and F 2 L , 

hong = -^=-V ub V* d (^C 2 + — bl J rUpE^'El' 



A lon g = -^V ub V: d (C 2 + —C 1 )m p eMeM 

{ - * [g»M ■ p) - p^u] ■ 2Fi(q 2 ) + wpV • 2F 2 L ( g 2 )} . (23) 



Again, one has to invoke a model to calculate the form factors. Estimates from the light-cone 
QCD sum rules give flT| : 



Ff/Fs = 0.0125 ± 0.0010 , F%/F s = 0.0155 ± 0.0010 , (24) 

where the errors correspond to the variation of the Borel parameter in the QCD sum rules. 
Including other possible uncertainties, one expects an accuracy of the ratios in ( f2~4"|) of order 
20%. Since the parity-conserving and parity-violating amplitudes turn out to be close to each 
other, F± ~ F 2 L = Fl, the ratio of the LD- and the SD- contributions reduces to a number 

A I A — r>B u ^P~f y ubv u d / 9 r\ 

s^longl •^-short — ^L/S \r \r* ' \ > 

VtbVtd 

Using C 2 = 1.10, Ci = -0.235, C e 7 ff = -0.306 (at the scale fj, = 5 GeV) § gives: 

B ^ A-n 2 m p (C 2 + CJN C ) F^ n . . 

R?^ n = pK „ ' ' ■ -4 = -0.30 ± 0.07 . 26 

L/S m b C e 7 ff F^^ 1 ; 

To get a ball-park estimate of the ratio Ai ong / A s hort, we take the central values of the CKM 
matrix elements, V ud = 0.9744 ± 0.0010 |2§, \V td \ = (1.0 ± 0.2) x 10~ 2 , \V cb \ = 0.039 ± 0.004 
and \V ub /V cb \ = 0.08 ± 0.02 0, yielding, 

\A long /A short \ B ^ = \Rfj^\\^M ~ 10% . (27) 

\VtdVbt\ 

The analogous LD-contributions to the neutral B decays B d — > P7 and B d — * w~f are expected 
to be much smaller, a point that has also been noted in the context of the VMD and quark 
model based estimates [3~1|. In the present approach, the corresponding form factors for 



the decays B d — > p°(u;)7 are obtained from the ones for the decay B u — > p ±/ y discussed 
above by the replacement of the light quark charges e u — > e d , which gives the factor —1/2; 
in addition, and more importantly, the LD-contribution to the neutral B decays is colour- 
suppressed, which reflects itself through the replacement of the factor a\ = C 2 + C\/N c in 
(p3|) by a 2 = C\ + C 2 /N c . This yields for the ratio 

= e J-l ~ -0,13 ± o,05 ; (28) 

where the numbers are based on using 02/01 = 0.27 ± 0.10 [[3^]. This would then yield at 
most R^ 1 - #f/<T 7 = °- 05 ' which in turn S ives A fon7 P yAsho7r < °- 02 - Even if this 



in 



underestimates the LD-contribution by a factor 2, due to the approximations made in [19, 20 
it is quite safe to neglect the LD-contribution in the neutral 5-meson radiative decays. 

The ratio of the CKM-suppressed and CKM-allowed decay rates for charged B mesons 
gets modified due to the LD contributions. Following [[3(J, we ignore the LD-contributions in 
T(B — > K*j). The ratio of the decay rates in question can therefore be written as: 

T(B U -* pry) 2 2 2i 



Ku \*[(l-p) 2 + v 2 ] 



T(B U -> K*j) 

x {l + 2 ■ RL/sV j^-A-f + {R L/ sYVt d {l P * + 2 ^ 2 | , (29) 

Using the central value from the estimates of the ratio of the form factors squared k u = 
0.59 d= 0.08 M, and the presently allowed range of the Wolfenstein parameters p and rj, it is 



shown in [19] that the effect of the LD-contributions is modest but not negligible, introducing 
an uncertainty comparable to the ~ 15% uncertainty in the overall normalization due to the 
S77(3)-breaking effects in the quantity k u . 

Neutral £>-meson radiative decays are less-prone to the LD-effects, as argued above, and 
hence one expects that to a good approximation the ratio of the decay rates for neutral B 
meson obtained in the approximation of SD-dominance remains valid 0: 

^A 2 [(l-p) 2 + ry 2 ] . (30) 



V(B -> K*-f) 



Here K d represents the SU(3)-breaking effects in the form factor ratio squared. It is a realistic 
hope that this relation is theoretically (almost) on the same footing in the standard model as 
the one for the ratio of the B°-B° mixing-induced mass differences, which satisfies the relation 



31 



AM, V u 2 



AM, 



Vt 



Id 



Ksd \2[(l-p)Z + V *} ■ (31) 



The hadronic uncertainty in this ratio is in the SU(3)-breaking factor K s d = (/l s ^B s //I d -B.B d ), 
which involves the pseudoscalar coupling constants and the so-called bag constants. This 
quantity is estimated as n sd = 1.35 ±0.25 in the QCD sum rules and lattice QCD approaches. 
(For details and references, see ||31|| ). The present upper limit for the mass- difference ratio 
AM s /AM d > 12.3 at 95 % C.L. from the ALEPH data |37|] provides better constraint on 
the CKM parameters, yielding |14d|/|14 s | < 0.35 than the corresponding constraints from the 
rare radiative decays B — > (p,uj) + 7, which give an upper limit of 0.75 for the same CKM- 
ratio. We expect experimental sensitivity to increase in both measurements, reaching the 



level predicted for this ratio in the standard model, | V^|/| V^,] = 0.24 ± 0.05 [|3T|, in the next 
several years in the ongoing experiments at CLEO, LEP and Tevatron, and the forthcoming 
ones at the B factories and HERA-B. 



Finally, combining the estimates for the LD- and SD-form factors in 19 
spectively, and restricting the Wolfenstein parameters in the range —0.4 < 
0.2 < 7? < 0.4, as suggested by the CKM-fits |l 
absolute branching ratios: 



and P, re- 
p < 0.4 and 
we give the following ranges for the 



B{B U -^ P1 ) = (1.9 ±1.6) x HT 6 , 

B(B d -> pi) ~ B(B d uj-i) = (0.85 ± 0.65) x 10~ 6 



(32) 



1 1 



where we have used the experimental value for the branching ratio BiB — > K* + 7) = 
(4.5 ± 1.5 ± 0.9) x 10~ 5 |Ij], adding the errors in quadrature. The large error reflects the 
poor knowledge of the CKM matrix elements and hence experimental determination of these 
branching ratios will put rather stringent constraints on the 

Summarizing the effect of the LD-contributions in radiative B decays, we note that they 
are dominantly given by the annihilation diagrams. QCD sum-rule-based estimates are very 
encouraging in that they lead to the conclusion that such contributions are modest in exclusive 
radiative B decays, in particular in the neutral 5-decays B° — > (p°,uj) + 7. This estimate 
should be checked in other theoretically sound frameworks. Of course, forthcoming data on 
specific i?-meson decays will be able to check this directly. Presently available data suggest 
that the contribution of annihilation diagrams in B decays is not significant, as seen through 
the near equality of the lifetimes for the B ± , B® and B® mesons and the near equality of the 
observed B^ and B° radiative decay rates. We have argued that this is very probably also 
the case for the CKM-suppressed radiative decays, with B ± — ► p ±/ y modified by O(20) % from 
its SD-rate. 
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